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We demonstrate the feasibility to completely characterize entanglement by negativities of 
quasiprobabilities. This requires the complete solution of a sophisticated mathematical problem, 
the so-called separability eigenvalue problem. Its solution is obtained for a non-Gaussian continu- 
ous variable quantum state, a two-mode squeezed state undergoing dephasing. This is a standard 
scenario for experiments with quantum-correlated radiation fields. 
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I. INTRODUCTION 

Entanglement is a key resource of quantum technology, 
for an introduction see [IHS] . The nonclassical properties 
of entangled quantum states have been discussed since 
the early days of quantum physics [U [5] . Applications 
of entangled states are, for example, quantum key dis- 
tribution [6], quantum dense coding [7], and quantum 
teleportation [5]. 

The most prominent and frequently used continuous 
variable (CV) entangled quantum state is the two-mode 
squeezed-vacuum state. It has been considered, for ex- 
ample, as a resource for quantum teleportation [51 ITU], 
quantum dense coding , and quantum memories |12j . 
It has been shown that this state violates a continuous 
variable Bell inequality [15] , Noise and loss effects of 
pure two-mode squeezed-vacuum states have been stud- 
ied, such as noise due to the transmission in optical 
fibers |14j . or dephasing and amplitude damping [15 . 

In general, nonclassical effects in Quantum Optics 
are characterized by the most prominent example of 
a quasiprobability distribution, the Glauber-Sudarshan 
P function [TUl [T7] . Negativities of this distribution ver- 
ify that the corresponding quantum state cannot be inter- 
preted as a classical mixture of coherent states, the latter 
being the closest analog to the classical behavior. Some- 
times the negativities are hidden in the highly singular 
structure of the P function. Recently, it has been shown 
in theory [T5] and in experiment [19 that these negativ- 
ities can be revealed by filtered, regularized quasiproba- 
bilities. 

The problem to verify the property of entanglement 
of a quantum state is very complex in general. Most 
difficult is the situation for non-Gaussian CV entangled 
states. One possibility is the use of matrices of higher- 
order moments [2D], which has been successfully applied 
recently [21]. Alternatively, entanglement witnesses can 
identify any kind of entanglement [221 E3]. However, 
the first method needs to identify particular matrices 
of moments out of an infinite manifold. In the second 
method one has to find the optimal witness among an 
uncountable number of Hcrmitian operators. Both types 
of entanglement tests are very difficult for non-Gaussian 
states. 

In the present article, we aim to provide methods of 



quasiprobability distributions for the property entangle- 
ment. Explicit solutions are derived for phase-diffused 
two-mode squeezed vacuum states, which are nontriv- 
ial examples of non-Gaussian CV entangled states. Our 
approach is based upon the previously developed con- 
cept of entanglement quasiprobabilities which requires 
to solve the so-called separability eigenvalue (SE) equa- 
tions for the quantum state under study [23] . Negativi- 
ties of these quasiprobabilities for a given state have been 
shown to be equivalent to entanglement. The general so- 
lution of this problem will be derived for our example 
of a class of non-Gaussian CV entangled states, which is 
of great interest for applications in quantum technology. 
Our method solely requires the quantum state and the 
solution of its SE problem, avoiding the complex iden- 
tification of moments or witness operators. Note that 
the connection between quasiprobabilities and their ap- 
plications to quantum technology has been reviewed in 
Ref. [23- 

In Sec. [Tl] we consider a realistic scenario for the gen- 
eration of the considered kinds of states. We recapitulate 
the basic reconstruction procedure for the entanglement 
quasiprobability in Sec. |III| The general solution of the 
separability eigenvalue equations is derived in Sec. |IV| 
On this basis we study in Sec. [V] the resulting entan- 
glement quasiprobabilities, which visualize entanglement 
via their negativities. A brief summary and conclusions 
are given in Sec. [VI] 



II. EXPERIMENTAL SITUATION 

In the following we consider the generation and the 
characterization of a nontrivial class of bipartite entan- 
gled CV states. They are created via phase random- 
ization of Gaussian two-mode squeezed-vacuum states, 
which eventually leads to so-called non-Gaussian states. 
They are not pure states anymore, but an infinite mixture 
of pure states. Altogether the states have complex prop- 
erties, so that an entanglement test is a demanding task. 
In the following we will not only witness the entangle- 
ment, but we visualize the entanglement by negativities 
of entanglement quasiprobabilities. The latter are nec- 
essary and sufficient for entanglement and contain the 
complete information on the quantum state. 
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Let us consider the generation of a two-mode squeezed- 
vacuum by an optical parametric amplifier |26) . 



£o = (i-c 2 ) E c 



m, m)(n, n|, 



(i) 



which is a pure entangled state and £ is the squeezing 
parameter (0 < £ < f). Let us consider an additional 
phase randomization in one of the channels [37], cf. the 
scheme in Fig. [I] We may model the phase randomization 
by a local unitary phase transformation 1 A £g) U(ip) which 
is performed with a certain probability p a (ip). As an 
example, we will consider a 27r-periodic Gaussian phase 
distribution with variance a, 



U(cp) = exp [irrup] 



= E 



exp 



2<T 2 



V2 



(2) 
(3) 



Due to the central limit theorem, such a distribution 
may often occur in practice, whenever small dephasing 
effects are added up. The resulting phase randomiza- 
tion probability is a Gaussian distribution with respect 
to the phase (p. In phase space (position-momentum rep- 
resentation), however, Po-(v) is not a Gaussian function 
anymore. 




FIG. 1: A two-mode squeezed- vacuum state is generated by 
an optical parametric amplifier (OPA). One mode is phase 
randomized by a random change of the optical path length. 
Correlation measurements are performed in both modes (Da 
and Db) for a reconstruction of the density matrix in Fock 
basis. 



Such a state Q a is, in general, experimentally accessible, 
cf. [26l[27] . For related states, the quantification of entan- 
glement has been investigated theoretically (35] HH] > and 
entanglement tests have been experimentally performed 
in reduced subspaces 2(g) 2, for detail see [3D] . However, in 
a 2 <gi 2 Hilbcrt space the Peres-Horodecki entanglement 
test is necessary and sufficient [33ll31| . In this paper, we 
go beyond this limitation. Moreover, the entanglement 
test is given in a form that includes the full information 
on the quantum state. 



III. DETERMINATION OF ENTANGLEMENT 
QUASIPROBABILITIES 

A first indication of the existence of an entanglement 
quasiprobability was the proof of the following represen- 
tation of general entangled states p in terms of separable 
ones: 



p = (1 + A)p scp - \p' s 



(5) 



scp ' 

with A > 0, p SC p and p' sep being separable states [32|l33]. 
Recently, we have shown how to construct entanglement 
quasiprobabilities on this basis [24 . This method en- 
ables us to show that an entangled quantum state p can 
be written in terms of separable states \a, b)(a, b\ and a 
particular quasiprobability PEnt (a, b), 



dPE nt {a, b)\a 7 b)(a,b\, 



(6) 



in analogy to the Glauber-Sudarshan P function. It is 
worth to mention that we can always rewrite the integral 
as an infinite (sometimes finite) sum, 



k 



PEnt{ak,b k )\a k ,b k )(a k ,b k \, 



(7) 



with some PEnt(afc, ^fc) < for an entangled state p. 

The initial difficulty in defining such an entangle- 
ment quasiprobability PEnt was the mathematical prob- 
lem that the decomposition of general states in terms 
separable ones is not unique. This means that a non- 
optimized decomposition of a separable state may con- 
tain negativities. Our method overcomes the ambiguity 
of the decomposition by optimization, based on the solu- 
tion of the separability eigenvalue problem of the quan- 
tum state. The state is entangled, if and only if the 
optimized PEnt contains negativities. 



The dephased squeezed state is given by 



Qc = 



d<p P(T (<p) {I A <8 U(cp))go(I A ®U(<p) f ) 



a 2 (m — n) 2 



(4) 



m,nSN 



x exp 



|m, m)(n, n\ 



A. Separability eigenvalue equations. 

Before we consider the reconstruction of PEnt, let us 
recall a mathematical tool which is needed to obtain the 
optimized entanglement quasiprobability. This tool is the 
set of separability eigenvalue (SE) equations [331 121] , 



Pb\a) = g\a) and p a \b)=g\b), 



(8) 
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with the reduced operators p a = Tr^([|a)(a| <S) and 
Pb = Tr b([^a ® |6)(6|]/o), with E^ ( b being the unity oper- 
ator in subsystem A,B. Using these equations one can, 
for example, construct all bipartite entanglement wit- 
nesses [23] . Applying the SE equations to a rank one 
operator [23], e.g. \tp}(tp\, yields the well-known Schmidt 
decomposition of l^) [3J. The value g in Eq. ^ is referred 
to as separability eigenvalue (SE value), and the vector 
\a, b) denotes the separability eigenvector (SE vector). 

The notion and properties of the SE equations and its 
solutions are in full analogy to the well-known ordinary 
eigenvalue equations. Solving the ordinary eigenvalue 
problem, p\ipk) — Pk\^Pk), enables us to find the spec- 
tral decomposition of a quantum state p in terms of pure 
states IV'/cKV'fcl and non-negative eigenvalues p k . When- 
ever an eigenvalue is negative, the corresponding oper- 
ator cannot be a valid density operator. In the present 
case we obtain the optimal representation of p in terms of 
product states defined in Eq. (JgJ) . Whenever the weight 
factor Pent contains negativities, the state cannot be sep- 
arable; hence it is entangled. 



B. Reconstruction of PEnt- 

As a fundamental application of the SE equations, let 
us study a three step reconstruction scheme to obtain 
fEnt, cf. [21] ■ In a first step, we solve the SE equation 
for the quantum state p, 

(i) Pb k \a k ) = g k \a k ) and pa h \h) = 9k\h)- 

Note that all SE values g k are non-negative, and a direct 
decomposition of p as in the case of the spectral decom- 
position is impossible. In a second and third step we 
can obtain the entanglement quasiprobability PEnt- We 
define 

(ii) the Matrix G = (\{a k , b k \a h h)\ 2 ) k . u 
a vector of the SE values g = (g k ) k , 

and a quasiprobability vector p = (pEnt(az.> 

Using these definitions, we simply invert the linear equa- 
tion 

(iii) Gp = g 

to obtain PEnt- Finally, the quantum state under study 
can be written as given in Eq. ([7]). 

Performing the steps (i) - (iii), we obtain the optimized 
PEnt of a given state p. At least one negative value in 
the determined PEnt is necessary and sufficient to verify 
that the corresponding quantum state is entangled. In 
analogy to the spectral decomposition, the quantum state 
under study is not a valid separable state. 



IV. SOLUTION OF SE EQUATIONS 

For general quantum states the solution of the SE 
equations is a cumbersome mathematical problem, for 



which no systematic methods are available yet. In this 
section, we solve this problem for a class of density oper- 
ators of the form 



P 



P m ,n\m,m)(n,n\, 



(9) 



with p m ^ n = p* nm - To this class belong the two-mode 
squeezed-vacuum states and their dephased versions. 
The solution is a substantial progress, since among the 
practically used CV quantum states the squeezed ones 
play a dominant role. The characterization of dephasing 
includes important effects which typically occur during 
the propagation of the radiation field. 



A. General Solution 

We seek the general solution of the SE equations for 
quantum states of the type given in Eq. The reduced 
operators for the states 



= J> fc |fc) and \b)=Y,h\l) 



(10) 



can be written as 
p a =Tr A ([\a)(a\ ®Ib)]p) = ^ a mPm,na n \m)(n\, (11) 

Pb =Tk fl (pr A ® \b)(b\)]p) =Y,KnPm,nb n \m){n\. (12) 

m,n 

The formulation of the SE equations ([8| in Fock basis 
yields, for p a \b) = g\b): 



and for p \a) = g\a) 



i)=^2gb m \m), (13) 



\m) =^2ga m \m). (14) 



Equivalently, for all m S N holds 



* T 

a m 9 ®rn ) 
m ~~ 9 ®"mi 



(15) 
(16) 



with normalizations X^gn l a «l 2 = SneN \K\ 2 = I- 

This system of equations can be solved by inserting the 
equations into each other. We obtain 



g 2 a m = 0, (17) 
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and analogous relations for b„ 
given. For all m £ N holds 



Now, the solution can be 



'm,n u "n u n 



9 or 



(18) 



which already delivers all solutions. This given form of 
the solution can be interpreted as a relation between dif- 
ferent coefficients of the solution. Either the projections 
|(m, m|p|a, b) \ = g are equal for different values of m, 
or the component does not contribute to g at all, for 

Q"m : bra — 0. 

In the following, we rewrite the general solution in a 
more convenient form for practical calculations. For this 
purpose we introduce the set Af which includes all indices 
m fulfilling the second part of condition (18), 



Af = {m € N 



b,, 



0}. 



(19) 



The remaining indices are elements of the set Af = N\Af. 



For those elements the first part of condition ( 18 1 can be 
written as 



neJ7 



'in ,n u ii u n 



(20) 



for arbitrary phases < <j) m < 2tt. In addition, we intro- 
duce some abbreviation for two vectors c and e, and the 
coefficient matrix pjf, 



e = (e"* m ) me AT, 



(21) 



PjT — {Pm,n) r 



Note that we have changed the normalization of the sepa- 
rable states, a m b m — > a m b m /g. This enables us to get rid 
of g. A trivial normalization of the unormalized vectors 
at the end of the procedure delivers the correct states. 



We simply rewrite Eq. ( 20 1 as 



(22) 



We obtain our SE vectors from the relation c m — a m b m , 
including the normalization. The SE values follow from 
the fact that g = (a,b\p\a,b), cf. Eq. We can 

choose now arbitrary subsets Af of all non-negative in- 
tegers N, phase vectors e and b m — c m /a m (normalized 
and m G Af). Each choice will deliver a particular solu- 
tion. This highly degenerate problem can be restricted 
to some principal solutions. 



B. Principal Solutions 

It is useful to restrict the manifold of solutions to the 
reasonable ones. For this reason we consider properties of 
our operator p. In particular, the operator p is symmetric 
(exchange of the subsystems A and B delivers the same 



operator), and it consists only of real variables, p m ,n = 
p* mn - Thus, we may assume \a, b) = \a, a), cf. a similar 
argumentation in [331 [3S], an d only real vectors e. This 
particular choice will be sufficient for the construction of 
the desired quasiprobability of entanglement. Using the 
solutions given above we can now choose a set Af, which 
yields 



i) = 

in 

= (a, a|L|a, a), 



1/2 



(23) 



e 



being the possible roots for each component, and 
= i e n) ne jf with e n <= {+!> This form is useful for 
numerical purposes, since it only requires the inversion 
of the N x N matrix pjj, where N = \Af\. 

Let us comment on the trivial solutions. Condi- 



tion ( 18 ) is fulfilled in the trivial case when for each m 
holds either a m = 0, or b m = 0, or a m = b m = 0. From 
a m bm = follows that the SE value is g = 0. Note that 
these solutions include all SE values g — 0, if we assume 
that the coefficient matrix (p m ,n)m,neN has full rank. 



C. Elementary Examples 

Example N = 1. For a better understanding of the 
solutions we may study the case Af = {k}. We insert this 
in the above general solution. This leads to a SE vector 
|a,a) = |fc, k) in Fock basis and a SE value g — pk.k- 
However, we can directly check whether this is a solution. 
The reduced operators used in the definition of the SE 
Eq. ([8) are 



Pa = Pb = Pk,k \ k)(k\. 



(24) 



Obviously, |fc) solves the eigenvalue problem for each re- 
duced operator. It is noteworthy that these solutions for 
Af = {k} already include the maximal SE value. This 
maximal SE value is needed to construct an optimized 
entanglement witness from the operator p, cf. [23l [28] . 

Example N — 2. A nontrivial example can be given 
for a set Af = {k, I}. This yields a matrix 



P{k,i} 



P{k\i} 



Pk,k Pi^ k 
pl,k Pi, I 
1 

dct (P{k,l}) 



(25) 



Pi,i 
-Pi.k 



~Pl,k 

Pk,k 



Using the definition of e = (efc,e;) T , the SE vector \a, a) 
and SE value g can be directly formulated as 



y/pue^-pik^lty ± yf-Pi,^" + Pk,k^ l \l) 



\pl,ie^" - pl k e^ | + | - pi.^fx + Pk^e 1 "" 
Pk,k + Pi.k^^ + Pike~ i{4 ' k " t " ) + Pl,l 



Kze* - p* lk e^ | + | - pi, k e^ + p k .k^ 1 \ ' 



(26) 
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with using the fact that a global phase can be ignored. 
Note that a global phase transformation can be per- 
formed for the state \a). The solution in this case only 
depends on the phase difference — 

V. APPLICATION TO DEPHASED TWO-MODE 
SQUEEZED- VACUUM STATES 

In this section we apply our method to the dephased 
two-mode squeezed vacuum state g a given in Eq. Q. 
For the identification of entanglement of a general CV 
entangled state it has been proven that it is necessary and 
sufficient to identify entanglement in finite subspaces of 



the Hilbert space of the system under study, for details 
see [36j . This is very helpful in practice. Making use 
of this result, entanglement of a non-Gaussian bipartite 
CV quantum state could be already demonstrated in two- 
qubit subspaces |30| . 

Based on this knowledge, we start to consider the en- 
tanglement quasiprobabilities in the Fock spaces contain- 
ing zero and one photon per mode. This subspace has the 
highest probability to be measured, and the interference 
terms are the maximal ones, cf. Eq. @. In Fig. [2] we il- 
lustrate the phase distributions chosen for our examples. 
They are ranging from a S distribution for a = towards 
an almost uniform phase distribution for a = 5. 




2 2 22'22 22 



FIG. 2: The distributions p&(4>) are shown for the considered phase randomization, cf. Eq. J3L for a — 0,0.5,2,5 (from left 
to right). The first plot represents a 8 distribution centered at the origin. The last plot is close to a full randomization with an 
almost uniform phase distribution. 



In Fig. [3] we show the entanglement quasiprobabilities 
for the two-qubit spaces under consideration. The first 
three distributions show clear negativities, which tend to 
decrease with increasing dephasing. These negativities 
uncover the entanglement of our dephased squeezed vac- 



uum state g a , not only for the chosen subspace but even 
for the complete CV state. For a full phase randomiza- 
tion the negativities disappear, since the corresponding 
state is separable. The plot for a = 5 is very close to this 
situation - the negativities are too tiny to be visible. 




FIG. 3: The plot shows the optimized entanglement quasiprobabilities for different values a = 0, 0.5, 2, 5 of dephasing, for 
a squeezed vacuum state with £ = 0.62, corresponding to —6.3 dB noise reduction. The ordering is the same as in Fig. [2] 
The values of k number the different solutions of the SE equations. The corresponding sets Mk are: Mi = {0}, M2 = {1}, 
A/"3,4,5,6 = {0, 1}. For the latter set the individual solutions are discriminated by different choices of the signs ±. 



r 



In Table [T] we give the SE values, the coefficients of the SE states, and the quasiprobabilities, for the case 



G 



of dephasing with (7 = 2. The numerical values are ob- 



tained from Eq. ( 23 1 . Taking into account that a global 



phase can be ignored, we obtain the given six SE vec- 
tors. Comparing the reconstructed state £> re c,<T according 
to Eq. ([7]) with the original, but truncated state g a as 



given in Eq. Q yields a numerical error of e = 4 x 1CP 16 
(precision of double-floating-point numbers). The nu- 
merical error is computed in the Hilbert-Schmidt norm 

e = [Tr(g Icc ^ - p CT ) 2 /Tr^] 1/2 . 



k 


1 2 3 4 5 6 


SE value gu 


0.615600 0.236637 0.190946 0.190946 0.149659 0.149659 


SE vector a^} 


1.000000 0.000000 0.496987 0.496987 0.549275 0.549275 
0.000000 1.000000 0.867758 -0.867758 0.835642i -0.835642i 




0.618287 0.217677 0.069431 0.069431 -0.061295 -0.061295 



TABLE I: The table contains the explicit numerical solution for a phase randomization a — 2. The columns include the 
information for the k-th solution. The two components ak,o and ak,i of the SE vector \ak) = afe,o|0) + afc,i|l) are given in the 
third and fourth line, respectively. 



In cases when the entanglement quasiprobabilities of 
the two-qubit subspaces do not show significant negativ- 
ities, we can proceed to study subspaces of higher dimen- 
sions. Let us consider two-qutrit subspaces including the 
photon numbers of up to two per radiation-field mode. 
In Fig. g] we deal with an example for relatively strong 
dephasing, a = 5. The entanglement quasiprobability 
in the two-qubit subspace fails to show significant neg- 
ativities, see Fig. [3j Our calculation yields a maximal 
negativity of —1.773 x 10 -6 . However, in the two-qutrit 
(3 (g> 3) subspace the entanglement quasiprobability be- 
comes clearly negative, with a maximum negativity of 
—3.786 x 10~ 2 . Surprisingly, even for such strong dephas- 
ing our method may identify some residual entanglement 
in the dephased two-mode squeezed vacuum state under 
study. In principle our method can be further extended 
to higher dimensional subspaces whenever needed to in- 
crease the significance of the negativities, at the expense 
of increasing numerical complexity. 




5 10 15 20 25 30 



FIG. 4: The quasiprobability of Q a is shown for a — 5 and 
a squeezing parameter £ = 0.62. The values of k number 
the different solutions of the SE equations, according to the 
sets: A7i = {0}, AT 2 = {1}, Jf 3 =_{2}, Ma,.. .,7 = {0,1}, 

ATs,...,!! = {0, 2}, ATl 2 ,...,15 = {1, 2}, ATl5,...,31 = {0, 1, 2}. 



Let us conclude this section with a general remark 
concerning the use of entanglement quasiprobabilities. 
We have outlined in Sec. |III| that the entanglement 
quasiprobabilities yield a complete representation of a 
general bipartite quantum state, cf. Eq. Q. In the suc- 
cessive treatment of subsystems of increasing dimensions, 
however, of course we only get a complete representation 
of the reduced quantum state in the chosen subspace. 
Another important issue is the quantification of entan- 
glement. In a very general sense entanglement can be 
quantified by the Schmidt number (SN), cf. [2^ 1 137 1 155], 
A method for determining the SN has been introduced 
recently [28] . In fact, the quasiprobability method devel- 
oped here can be extended to identify the SN via nega- 
tivities of SN-quasiprobabilities. 



VI. SUMMARY AND CONCLUSIONS 

We have derived the entanglement quasiprobabilities 
for a nontrivial entangled quantum state. This distri- 
bution visualizes nonclassical correlations between radi- 
ation modes. Our studied state consists of a two- mode 
squeezed vacuum undergoing dephasing in one channel. 
The resulting state is non-Gaussian, which is of some 
interest for various applications in quantum technology. 
We provide a method for solving the separability eigen- 
value problem of the quantum state under study. From 
the solution we could derive the entanglement quasiprob- 
abilities of the state. Their negativities are necessary 
and sufficient for entanglement of the truncated quan- 
tum state in a chosen subspace. Moreover, the existence 
of negativities in some subspace is necessary and suf- 
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ficient for the entanglement of the infinite dimensional 
continuous- variable quantum state. 

We have identified entanglement of dephased squeezed 
states in two-qubit subspaces. For strong dephasing, 
however, the corresponding entanglement quasiprobabil- 
ities may only show tiny negativities which are hard to 
determine in experiments. In such a case we have shown 
that the entanglement quasiprobabilities in a two-qutrit 
subspace may show negativities which are enhanced by 
a factor of about 2 x 10 4 compared with the two-qubit 
space. Beside the more significant entanglement signa- 
tures, the extension of our method to higher dimensions 
is of fundamental interest. For this purpose, we formu- 



lated the general procedure and the solution of the sep- 
arability eigenvalue equations. It allows one to identify 
quantum correlations for the important class of radiation- 
field states under study, in the notion of bipartite entan- 
glement in its most general form. 
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